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Consider a deformable body immersed in an incompressible fluid that is 
randomly stirred. Sticking to physical situations in which the body departs 
only slightly from its spherical shape, we investigate the deformations of the 
body. The shape is decomposed into spherical harmonic modes. We study the 
correlations of these modes for a general class of random flows that include 
the flow due to thermal agitation. Our results are general in the sense that 
they are applicable to any body that is described solely by the shape of its 
surface. 

I. INTRODUCTION 

The current paper is the second in a series of papers that treat systems of deformable 
objects immersed in an incompressible liquid. These systems are common in nature ( milk, 
emulsions, composite fluids )and have been studied extensively. Nevertheless, there are 
many questions concerning them that need to be answered still. This is so, because the 
actual solution of such systems is extremely difficult. The objects interact with themselves 
and with each other via the hydrodynamic interactions, and form, in general, a non-linear 
many-body problem. Moreover, each object has infinite degrees of freedom that correspond 
to deformations, and all of them , of all the objects, interact. The problem is a little simpli- 
fied if the deformation of the objects from spherical shape remains small and the fluid is in 
the linear regime, as the case of the Stokes approximation to the Navier - Stokes equation. 
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Our final goal is to obtain the response of the composite system to a given velocity field 
imposed on the liquid. The velocity field we have in mind may be fixed in time like sim- 
ple shear or randomly fluctuating in time and space. Even in the first case the velocity 
field, each object experience must have a random part due to the random passage of other 
objects nearby. Therefore, effects of randomness are important to the understanding of 
such systems. We wish to know first, what is the response of a single deformable object to 
the random velocity field, created by the other objects and the external source. Once the 
response of a single object is known, we can obtain the response of the full, many body, 
system by using the response of each body as a source of an additional velocity field. In 
the former paper [Jj we studied the motion of the center of a deformable object and derived 
the Mean Squared Displacement of its center. In this paper we investigate the deformation 
degrees of freedom, within the scope of the simplifications stated above, thus completing the 
description of the response of a single object to a random velocity field. Decomposing the 
deformation into spherical harmonic modes, we consider the correlations between deforma- 
tion modes. We find, among other things that different modes are decoupled and that the 
correlation function does not depend on the parameter m of the Y^ m spherical harmonics 
(following from spherical symmetry). We also obtain a method to calculate these correla- 
tions as a function of time, build typical drop shapes from the correlation functions, and 
discuss several interesting cases as thermal agitation. 

The plan of this paper is as follows: In section || we describe the system we have in mind 
and formulate the basic equations. In section |T| we introduce the deformation coefficients 
and obtain their correlations for the simple case where the external velocity is uncorrelated 
in time. General correlation functions are considered in section |IV|. We obtain the defor- 
mation correlations and make simplifications for equal times. In section |V| we discuss some 
properties of the correlation function of the external velocity for several cases, among them 
is the special case of thermal agitation. An algorithm for numerical computation of the 
correlations in the general case is described in the Appendix. 
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II. THE SYSTEM 



Consider a single deformable body immersed in a host fluid. 

1. The deformable body is fluid, in the sense that the velocity fields are well defined 
everywhere (both inside and outside the body). In particular, each surface element 
moves with the velocity of the flow at its position: 



3. The body is characterized by an energy that depends on its shape ( Changing the 
orientation or switching places of two surface particles while keeping the shape constant 
does not cause any energy). The shape of minimum energy is a sphere. The energy 
may be surface tension ||, Helfrich bending energy |7|,|8|], etc. Deformation of the 
shape change the energy, exerts force density in the liquid and therefore creates an 
additional velocity field, denoted v^. 

4. We investigate the regime where the hydrodynamic equations are linear in the velocity 
( i.e. a velocity field induced by several sources is equal to the linear sum of the velocity 
fields that each source induces separately). For instance, if the flow is governed by 
the Navier Stokes equation, then our assumption implies that the Reynolds number 
is small and that the stokes approximation is applicable. The linearity implies that 
in our system the actual velocity field is the linear sum of the imposed velocity field , 
v ex t ( where v ex t is the velocity field that would have been the actual velocity field if 
the body was absent), and the velocity field induced by the deformations, v^, 




(1) 



2. Both the body and the host fluid are incompressible, V • v — 0. 



V = Vext + *V 



(2) 



5. We consider cases in which the external velocity field is random, with an average of 
zero. The velocity correlation function is known and depends only on distance and 
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time difference. We also assume that the the external velocity is weak enough to allow 
the body to remain almost spherical 



III. WHITE NOISE FLOW 

Consider a spherical body which is slightly deformed. The equation 

iL + /(n,t)-i = o (3) 

defines its surface , providing for each spatial direction, fl, the distance, p = \ f— r |, of the 
surface from the center of the body, r . R is the radius of the undeformed sphere. f(fl, t) that 
parameterize the shape and is named the deformation function. The deformation function 
is decomposed into spherical harmonics, f(£l,t) = Ylil-oc, ^2 m =-l flm(t)Yi m (Q), where the 
deformation coefficients are denoted fi m {t)- Our goal is to obtain the correlations between 
deformation coefficients. 

The center of the body, f*o, is chosen to be the point around which the deformation coefficients 
with I = 1 vanish: / l m = 0. A different definition of the center will introduce three 
additional equations for the deformation coefficients with 1 = 1. We are not interested in 
those since in the first order the spherical harmonics with / = 1 describe a translation of 
the body. In previous papers we have used the same definitions to describe the motion 
of the center. Therefore, their results are consistent with the derivation done here and can 
be used. Let ip(r) be a scalar field, defined everywhere in such a way that the equation 
ip{r) = describes the surface of the body ||f§]. Straight forward manipulation of eq. ([[]) 
give rise to a continuity equation for ip, presented here in a coordinate system that moves 
with the center of the body: 

ip + • Vip = -{v ext - f ) • Vij) (4) 

A good candidate for ip is found in equation (|J), ijj = ^ + f(Q,t) — 1. Assuming that \v ext — r\ 
is small, the right hand side of eq. is, in the first order, equal to Q = j/5 • [v ext — ^]| 
(see ||), where p is a unit vector directed outwards from the center of the body. Since 



the minimum energy of the body is obtained for a spherical shape, the velocity induced 
by the body is zero when the sphere is undeformed. Therefore, the velocity tty must be, 
in general, a linear functional of the deformation f(Q,t) ( with corrections of higher order 
that we neglect). The term ■ \7ip is obtained in the leading order by taking on the 
original sphere and ty, to first order in the deformation. We are interested in the first order 
evolution of the surface, and the terms above are already linear in the deformation; Thus, 
we can calculate them on the original sphere, where it is possible to develop eq. (f|) in terms 
of spherical harmonics (for further explanation see |]]). Consequently, the generic eq. for 
the deformation coefficient, fi >m , must be of the form 

^^ + \ l fUt) = -Qi m (t), (5) 

where Q\ m is given by 
1 ( 

Qlm = ^jdn{p- Vext-r Y^Q)} (6) 

and v ext is evaluated on the surface of the body (and can be approximated as on the unde- 
formed body) at the direction of the spatial angle Q. 

It is convenient to write the correlation function of the external velocity field in the 
momentum space. This is so because the random velocity field is transversal when the fluid 
is incompressible. Consequently, in real space, the flow must always be correlated in a very 
complex way. On the other hand, in the momentum space we can simply use a projection 
operator on the transversal part of a general field: 



VextM ) = E - ^rj 



Uj(q), (7) 



where u is a general vector field and the bracketed term is the projection operator that 
removes the longitudinal part of u, and therefore yields a general transversal velocity field 

Vext- 

Next, the correlations of the external velocity are easily expressed using the correlations of 
the general field u. We are interested in cases were the system is isotropic, homogeneous 
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and stationary. In these cases, the general field must obey: 

(ui(q, t^Umip, t 2 )) = 5 lm 5(q + p)(f){q, t 2 - h), (8) 

where 5i m is the Kronecker delta, 5() is the Dirac delta function and </>() is a general function 
of q and At. In addition we have assumed that 

(«,(?,*)) = 0. (9) 

In the rest of this section we consider a frequently used family of random flows in which the 
external velocity is uncorrelated in time. 

<P{qM-t 1 )=^{q)5{t 2 -t l ) (10) 

Equation ([TTJ) is a reasonable representative of systems in which the smallest time scale is 
that of the random velocity. In these systems we can replace the exact correlation details 
by the effective delta function by integrating on time. 

Using Fourier transform and the definition of the correlations we can calculate the correlation 
of the velocity at two places on the drop, that are characterized by the directions f and r' . 
The calculation yields 

QiQj 



(vUr,h)vUr',t 2 )) = / d?e*W 



m)s(t2-ti), (ii) 

where v ex t{r,t) is the velocity at time t at place r on the surface and r = r (t) + Rf. The 
average and correlations of Q^ m follow easily from the previous equations. The average of 
Qi m is zero, 

(Qim) = 0. (12) 

The term p ■ r , in eq. (0), does not contribute to any component of Qi m except for those 
with 1 = 1. In addition, the center, r , have been chosen to be the point around which the 
three deformation coefficients, fi m , with I = 1 are equal to zero. Therefore Q\ m is zero, and 
for all other choices r can be dropped out of Qi m . 

The external velocity on the surface is approximated, in the first order of deformation, 



as the external velocity on the undeformed sphere Strait forward calculation of Qz,m 
correlations using its definition, eq. (|J), and the velocity correlations, eqs. (0-|lO|), yields: 



(Qlm(tl)Ql'm' (h)) — Qlml'm'S(t 2 ~ ^l) > 



(13) 



where 



Si 



<Mj_ 
Q 2 



(14) 



Finally, the correlations of the deformation coefficients, fi m , are obtained, using Eq. 
and the correlations of Qi m , by direct integration 

rt+At 



(flm(t)fl'm'(t + At)) t ^ x> = lim / dti / dt 2 

t^°°Jo Jo 



\QlmV>l)Ql'm' (h)) e 



Aj(ti-t)+A;/(t 2 -t-At) 



(15) 



where the limit t — > oo is taken to avoid the initial conditions that are imposed on a specific 
realization. 

The last equation combined with Eqs. ( |i~4|) and (|i~3|) yields 

e -A,/|At| 



(flm(t)fl>m>(t + At)) 



t^oo 



Q 



Iml'm' 



(16) 



A/ + Ay 

Since the system is invariant to time reversal it is obvious that there is no preference of 
over A/ in the exponent, and therefore the result must equal zero for I ^ I'. In fact, only for 
the terms in which I = V and m = —m' {Qi m (0)Qi' m '{0)) ^ 0. This follows from the fact that 
in the first order Q\ m depends only on the component of the external velocity that is normal 
to the surface. Therefore, its correlation function for two places at the same time must 
depend only on the angle between their direction, and can be expanded as YLi AiPi(cos(9)) 
where Pi is the Legendre polynomial. This can be turned into a sum of spherical harmonics 
using the partial waves expansion 
Arc 



Pi(cos(9)) = — - £ {~l) m Y lm {m -m(tf) 



(17) 



m=—l 



Now, it is obvious that the correlations of the deformation function, f(£i), must have the 
same form 

47T - 
l m=—l 

Furthermore, 

(fimfiw) = J if(n)f(rt)) Y^Y^in') cin dsn. (19) 

Therefore, the expansion is composed of terms in which I = I' and m = —m' only, as stated 
above. 



IV. GENERAL NOISE 

In many cases white noise correlations are not sufficient to describe what really happens 
in the liquid. Especially if the correlation time is of the order of other time parameters. Such 
is the case of a system of many droplets inside a host fluid. The random flow, a droplet is 
subjected to, results from the random motion and deformation of other droplets that pass 
nearby. It is obvious that in this case the approximation of the flow to be uncorrelated in 
time is not justified. Hence, we need to generalize the theory. 

Suppose that <p(q, At) is a general function of q and the time differences At. The corre- 
lations of the external velocity are now extended in time. In order to calculate averages on 
the droplet in different times we must now consider also the movement of the droplet. The 
definition of Qi ym © implies that the correlations of the normal component of the external 
velocity field on the surface of the body are 

{Qi m {ti)Qvm>{t 2 )) = -L / dQ f dQ'Y; m (Q)Y; m! (Q') £ [ 

i,j=x,y,z 

^(vUr^vUr'M)]. (20) 
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The correlation of the velocity at two points on the droplet, located in the directions f and 
f and measured at different moments, obviously depends on the displacement of the center, 
Ar*o. We calculate first the correlation for a general displacement of the center and then, 
average the result according to the probability of finding the center at each point. To do this, 
we first express the velocity correlations in means of the Fourier transform of the velocity, 
use equations (0),(H) and obtain 

(vUr,ti)vUr',h)) = J P(Afo)d(Af ) J d 3 q e -^ +R ^[5 l3 - ^|]0(g, At),(21) 

where Ar*o = fo(ti) — ^0(^2), At = £2 — t\, P(Aro) is the probability that the center will be 
displaced by Ar in the period of At and the summation over d(Af ) is taken on all the 
possible configurations. 

It is obvious now that averaging over the center displacements will effect only the term 
— iq ■ Af in the exponent since this is the only term that depends on Ar . Consequently, 

{vUrM)vUr'M)) = f d 3 q(e- i ^)e- i ^-^ R (5 ij - )0(g, At), (22) 

J q 

Assuming Gaussian distribution of the displacements of the center, 

In a previous paper [0], we considered the Mean Squared Displacement (MSD) of the center 
of a deformable body in a flow that is correlated in a general way. We found that the MSD 
in a period of time At is given by: 

/•At poo 2 

F(At) = 16n dt' q 2 dqe-- Fit ^{q,t'){ ]0 {qR)+j 2 {qR)f{At~t'), (24) 
Jo Jo 

where F(At) = ((Af*o) 2 ). Therefore the correlation of the external velocity at two points on 
the surface, characterized by the directions f and f', measured at two different times with 
time gap of At, is given by 

{vl xt {rM)vUrM)) = [ ^-^e"*^ - )0(g, At), (25) 

J u 
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Equations fl24|) and fl25|) enable us to calculate the correlations of Qi tTn (eq. (p0[)). The 



correlations of the deformation coefficients can be obtained from their basic equation (|5|) 
using the correlations of Qi, m , 

(fi, m (t)fvMt + At)) t _ = J q dt J q dt exp{Xlt + Xl/{t + At)) • M 

Equations (PU[), fl2~j|), (f23| ) and (f2B| ) form the calculation method for the correlations of the 
deformation coefficients. The method presented here may be hard to implement to numeri- 
cal calculations because of the many dimensions integration. In the appendix we describe an 
algorithm that reduces the integrations to be one dimensional, thus making the computation 
task easier. This algorithm was used to obtain the results presented in the following section. 
At same times, At = 0, the correlations above play an important role in the shape distribu- 
tion of the drop. They hold all the information needed to determine the distribution of the 
// jm 's. As was explained in the previous section, same time correlation (STC) are obviously 
non-zero only for / = I' and m = —m! . In addition we consider cases where (fi >m ) = 0. Thus 
(fi,m(t)fi-m(t)) give the variance of the deformation coefficient fi >m) and since fi >m with 
different /, m are uncorrelated we can use this variance to create typical shapes of drops 
(Fig. §. 

Equation (f26|) can be simplify a little in the case of STC 



(AmW/r^WU = r dt '^(0)QvMt')) eXPi A ^ + f p( Xl,t '\ (27) 

JO Al + Al' 



V. THE CORRELATION FUNCTION 

We wanted the method presented above to be applicable to different random external 
flows. Therefore, we have built the calculation method in such a way that the correlation 
function <p(q, t) must be given by hand. In some cases <p(q, t) can be calculated theoretically. 
In others, to obtain <f)(q,t) for a given system one can make an experiment. Exclude the 
deformable body from the system; measure the correlations in the velocity in different points 
at different times and deduce the correlation function <ft(q,t). To confirm the results we 
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present here just insert back the deformable body, measure the correlations and compare 
with the theoretical prediction. We can also try to use our method to obtain the correlation 
function itself. Of particular interest is the response of a system of deformable bodies to an 
external flow. The deformation of each body will induce an additional velocity field that will 
influence the others. Hence, the correlations in the velocity, each body is subjected to, can 
be treated in a self consistent manner, using the deformations as a source to the external 
velocity field as well as the field's product. 

In the rest of this section we will analyze examples of random flows. 
Consider first the special case of thermal agitation. In a previous paper |l] we have shown 
that the correlation function for the external velocity due to thermal agitation have the form 
4>(q,t) = (^)3^ ~t-; where rj is the viscosity of the fluid (notice that the correlations are long 
ranged due to the incompressibility of the fluid). Using this form with eq. ([16]) the STC of 
the deformation coefficient, f\ m , where / > 1 is given by 

KgT 1 

(fl,m(t) fl-m{t)) '^oo = ^>2^ Q _|_ 2){l - 1) 

Schwartz and Edwards H considered the special case of deformable a body in equilibrium 
at temperature T, using the Kirkwood equation. They found identical correlations in the 
deformations (using X^ m = Rfi, m ) ■ Their derivation, however, has been tailored for thermal 
agitation and cannot be expanded to take into account any other correlations in the host 
fluid. 

The shape of the deformable body is described by r = R(l + /(O)). Hence, eq. (p8|) 
implies that the magnitude of the deformation, Rfi <m , does not depend on the size of the 
droplet, i.e. the deformations should be as of an infinite lamellar surface. 



Thermal agitation belongs to a class of systems in which the decay time, which is the 
time that takes the external velocity field to loose memory and become uncorrelated, is the 
shortest time scale (Although other time scales T\ = y tend to zero for I — > 00 we can 
use this form to any finite set of spherical harmonics. In addition we must remember that 
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/ is limited by the fact that features on the sphere must be larger than the inter-atomic 
distance for the continuum description to hold). In these cases we can approximate the 
actual velocity field to be totally uncorrelated in time. 

0(g,t) = C0(g)5(t). (29) 

There are two possible scenarios: either 0(g) decay with a length scale, £ (at least one), 
or it is long-ranged. We limit our discussion to one power and one length scale and the 
generalization is trivial. In the first case 0(g) = g _a 0(£g), where has a cut off at £g = 1. 
There are two dimensionless parameters for the STC, ji\ — j and ^ — rs-<*\ 1 ■ m the case 
of a body characterized by surface tension energy, A/ oc ^ where A is the surface tension 
and rj is the viscosity, 

(AmW/Cm'^)^ = Fl,m,l',m' i^ R A- a y " ( 30 ) 

Fig. [3] and £| depict the STC dependence on / for two correlation functions and different 
parameters values where we changed /ii and /i2 by keeping £ = 1 and — 1 and changing 
R. As can be seen, there are two possibilities: either the STC with I = 2 dominates the 
curve or there is a maximum at Zo ~ j- It is easy to see that the decay, for I >> Zo is 
exponential. This suggests that there is a cutoff on the deformation coefficients at which the 
expansion can be terminated and therefore that the expansion we use here will be useful for 
systems in which /i! is small. 

Fig |l] depicts various shapes of a body characterized by surface tension under a random flow, 
given by 0(g, t) = C5(£g — l)8(t). As can be seen, the surface of the body develops bumps. 
It is obvious that the size of these surface features depends on the ratio /ii = (Fig. 
[I]). As 7 increases, different surface elements become less and less correlated. Therefore, we 
expect to see features of smaller size (which correspond, clearly, to spherical harmonics of 
higher order). The smallest features correspond to deformation coefficients with / ~ fii (or 
I = 2 if n x < 2). 

In the second case the parameter /ii can be dropped out and we are left only with /i 2 , as in 
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the case of thermal agitation. 

The correlation of the deformation coefficients for different times (fi t m(t)fi',m'{t + Ai)}f_» 
depend in addition on the dimensional parameter fi^ = XiAt exponentially as written in eq. 

m 







FIG. 1. A typical realization of a deformable body subjected to a random flow of the form 
= C8(q£ - l)5(t) with (A) R = 2£, (B) R = 4£, (C) R = 6£ and (D) R = 8f . 
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FIG. 2. The variance of the deformation coefficients, fi m , as a function of Z, for a typical decaying 
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For R » £ the curves coincide with the curves for thermal agitation. (Go = M2) 



Other classes of random flows include separable correlation functions in which there 
is a decay length scale, £, and decay time scale r: <j)(q,t) = Cq~ a <fi(£q)T(t/T), where 
(/j,m(t)/i' ) m'(*)) t ^oo = F i,m,i',m> f 5 A,r) , and flows of the form <j>{q,t) = Cq^^Tqt^) 

for which analog dimensional parameters can be written. 
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VI. SUMMARY 



We have built a method to calculate the correlations of the deformation coefficients 
that correspond to the decomposition of the shape into spherical harmonics. We did it in 
two stages. The first, for external velocity fields that are uncorrelated in time: eqs. ( p~4|) 
, ([TT)|). The second, for a general external velocity field that is correlated both in space 



and time: eqs. (f2U|), (pip, (PS] ) and . The algorithm was based on the assumption that 
the correlation of the velocity field in an identical system without the body are known. We 
discussed these correlations, used the results to build in a rigorous way typical surface shapes 
and considered the special case of thermal agitation. In addition, we pointed out that from 
our numerical results there there is a parameter fii = where deformation coefficients 
with / > fix seem to decay exponentially and therefore are essentially unimportant. This 
suggests that working with spherical harmonics to investigate systems of deformable bodies 
is extremely useful for cases where Hi is small. 

APPENDIX A: ALGORITHM FOR THE CALCULATION OF THE 
DEFORMATION CORRELATIONS 



The velocity correlations involve a three dimensional integration, and Eq. (H) uses an 
additional four dimensional integration. As we can see, the method presented above is very 
hard to use. Therefore we must find a way to lower the dimensionality of the integrals, 
the following algorithm illustrates a method to do so using the addition theorem (or partial 
waves expansion) . The result is a finite expression which is composed of sum of terms where 
each term involves only one dimensional integration. Unfortunately, although finite, this 
sum is too long to be presented here. 
The partial waves expansion is given by 

oo I 

£ £ (-0^,(^)1^(0^(0), (Ai) 

1=0 m=-l 



-iq-{Rp) 
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where Q q is the solid angle in the q direction and ji is the spherical Bessel function. Using 
the partial waves expansion with Eq. (|TT| ) yields 



1 m 1 i 2 m 2 ijYi imi (Q)Y l2Tri2 (Q ) , ( A2) 

lx,h,mx,m2 

where 

A hmihm2ij = (-i) h i h (4tt) 2 | d?q e~i F ^ j h (qR)j l2 (qR) ■ 



q 2 



0(M,At) (A3) 



Each coefficient Ai m ii m iij divides easily into two parts, one deals with the angular integration 
and the other with the radial integration. The angular part is easily integrated, by the use 
of the Clebsch-Gordan coefficients, to give an algebraic expression. Therefore to evaluate 
Mmi'm'ij one needs to calculate the radial integration which is one dimensional. 
We use Eq. ( |A2|) and obtain a new equation for the Qi m correlations. 



1217121 j 

i^jnown)] (A4) 

I 1 ,m 1 ,l2,m2 

The integrations over Q and Q' can be easily calculated once we take into consideration that 
unit vectors can be decomposed into spherical harmonics with / = 0, 1. The integration 
yields a finite sum of terms. Each term is composed by a known number, (, multiplied by 
the corresponding Ai ltmi j 2tm2ti j that involves only one dimensional integration. Hence, 

(Qi m (0)Qi W (At)) = J2 

h,l2,m 1 ,m2,i,j 

Cw 2 ,™«(MVX) / q'dqe-^^MqR^iqRWlql, At) (A5) 



In this way we formulated a method for the calculation of the correlations. The explicit 
expression is extremely long and therefore is not presented in this paper. 
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